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Abstract. We present an analysis of a 10-day continuous ASCA observation of 
the narrow-line Seyfert 1 galaxy IRAS 13224-3809. The total band (0.7-10 keV) 
light curve binned with 500 s reveals trough-to-peak variation by a factor > 37. 
Rapid X-ray variability with a doubling timescale of 500 s has also been detected. 
The soft (0.7 — 1.3 keV) and hard (1.3 — 10 keV) X-ray band light curves binned 
to 5000 s reveal trough-to-peak variations by a factor > 25 and ~ 20, respectively. 
The light curves in the soft and hard bands are strongly correlated without any 
significant delay. However, this correlation is not entirely due to changes in the 
power-law flux alone but also due to changes in the soft X-ray hump emission 
above the power law. The variability amplitude changes across the observation 
but is not correlated with the X-ray flux. The presence of a soft X-ray hump 
below ~ 2 keV, previously detected in ROSAT and ASCA data, is confirmed. 
Time resolved spectroscopy using daily sampling reveals changes in the power-law 
slope, with Tx in the range 1.74 — 2.47, however, day-to-day variations in Tx 
are not significant. The Soft hump emission is found to dominate the observed 
variability on a timescale of ~ a week, but on shorter timescales (~ 20000 s ) the 
power-law component appears to dominate the observed variability. Flux resolved 
spectroscopy reveals that at high flux levels the power law becomes steeper and 
the soft hump more pronounced. This result is further confirmed using an earlier 
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ASCA observation in 1994. The steepening of the photon-index with the fluxes in 
the soft and hard bands can be understood in the framework of disk/corona models 
in which accretion disk is heated by viscous dissipation as well as by reprocessing 
of hard X-rays following an X-ray flare resulting from coronal dissipation through 
magnetic reconnection events. Time dependent accretion disk-corona models are 
required to understand the observed correlation between the soft hump emission 
and the power-law flux. 

Key words, galaxies: active - galaxies: individual (IRAS 13224-3809) - galaxies: 
nuclei - galaxies: Seyfert - X-rays: galaxies 

1. Introduction 

Seyfert 1 galaxies are an important class of active galactic nuclei (AGN). They show a 
large range in the width of their optical emission lines e.g., full width at half maximum 
( FWHM) of the H/3 line is found to be in the range - 1000 - 10000 km s" 1 . Seyfert 1 
galaxies that are at the lower end of the line width distribution with FWHM (H/3) Js 
2000 km s _1 are called the narrow-line Seyfert 1 galaxies (NLS1; Osterbrock & Pogge 
1985; Goodrich 1989) and are distinguished from the bulk of the Seyfert Is ("broad-line 
Seyfert Is" or BLSls). The NLS1 galaxies are also found to have strong Fe II emission and 
[O III]A5007/H/3 < 3 (Osterbrock & Pogge 1985; Goodrich 1989). However, Rodriguez- 
Ardilla et al. (2000) have shown that the ratio [O III]A5007/H/3 does not clearly distin- 
guish between NLS1 and BLS1. Also Veron-Cetty, Veron & Goncalves (2001) have shown 
that when only the narrow component of H/3 is considered, the above ratio is similar be- 
tween NLSls and Seyfert 2 galaxies. X-ray properties of NLSls are even more remarkable. 
These AGNs very frequently exhibit rapid and/or large amplitude variability (Boiler et 
al. 1996; Forster & Halpern 1996; Molthagen et al. 1998). The excess variance for NLSls 
is typically an order of magnitude higher than that observed for samples of BLSls with 
similar luminosity distribution (Leighly 1999a; Turner et al. 1999b). Giant-amplitude 
X-ray variability (up to a factor of 100) has also been observed in several NLSls (Boiler 
et al. 1997; Brandt et al. 1999). Some NLS1 galaxies show extremely rapid variability 
(on timescales of a few hundred seconds) by a factor of about 2-3 (Rcmillard et al. 1991; 
Boiler et al. 1997; Brandt et al. 1999; Dewangan et al. 2001a). ROSAT (0.1 - 2.4 keV) 
observations have revealed that the soft X-ray continuum slopes of NLSls are system- 
atically steeper than those of BLSls (Boiler et al. 1996), the photon index T x (photon 
flux f E oc E~ Tx ) sometimes exceeding 3. ASCA observations have shown that the hard 
X-ray 2 — 10 keV continuum slope too is significantly steeper in NLSls than that in the 
BLSls (Brandt et al. 1997; Turner et al. 1998; Leighly 1999b; Vaughan et al. 1999b). 
The very strong anti-correlation between FWHM of the H/3 line and both the X-ray 
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slopes in Seyfert Is (Boiler ct al. 1996) and in quasars (Laor et al. 1997), and "excess 
variance" (Turner et al. 1999b) suggests that the remarkable X-ray properties of NLSls 
are possibly due to an extreme value of a fundamental physical parameter related to the 
accretion process. 

A popular explanation for the distinct properties of NLS1 galaxies is that they have 
lower black-hole masses than the BLS1 galaxies. Smaller black-hole masses result in 
shorter timescales, thus naturally explaining the rapid X-ray variability, since the primary 
emission would originate in a smaller region around the central black-hole. Smaller black- 
hole masses also naturally result in narrower optical emission lines provided the size of 
the broad emission line region (BLR) scales with the luminosity (Laor 1998). 

A comparison of the soft X-ray properties of Seyferts and Galactic black-hole can- 
didates (GBHCs) led Pounds et al. (1995) to make an analogy between the two types 
of objects. They suggested that NLSls are the high state analog of BLS1 galaxies. The 
high state GBHCs show strong soft X-ray excess, with blackbody temperature ~ 2 keV, 
above a steep power law and are thought to emit a higher fraction of their Eddington 
luminosity. This led Pounds et al. (1995) to postulate that NLSls must also be emitting 
a higher fraction of their Eddington luminosity, hence higher accretion rates relative to 
the Eddington accretion rate (to = ^ M ) are required. Since NLSls have comparable 
luminosity to that of BLSls, a higher fractional rate also means a lower black-hole mass. 

The higher the fraction of the Eddington luminosity emitted, i.e. the higher the 
fractional accretion rate, the greater the temperature attained by the accretion disk, i.e. 
the disk emission becomes energetically dominant in the soft X-rays (Ross et al. 1992). 
Thus NLSls might be expected to show disk components which peak at higher energies 
than for BLSls. The spectral energy distribution (SED) from far-infrared to X-rays of 
NLS1 galaxies appears to be similar to that of BLS1 galaxies, but the UV luminosity 
of NLSls tends to be smaller than that of BLSls (Rodriguez-Pascual et al. 1997). The 
lower UV luminosity of NLS1 galaxies compared to that of BLSls could be due to the 
shift of the big blue bump (BBB) towards higher energies. The steep soft X-ray spectrum 
could be the high energy tail of the BBB (Mathur 2000). Pounds ct al. (1995) noted that 
the excess soft X-ray emission of NLSls may cause an increased Compton cooling of hot 
electrons in the corona resulting in a steeper hard X-ray power law. Higher accretion 
rates also result in an ionized surface for the accretion disk (Matt et al. 1993). Evidence 
for the ionized disk is found in the form of Ka emission from the ionized states of Fe 
in NLSls (Comastri et al. 1998; Turner et al. 1998; Vaughan et al. 1999a; Comastri et 
al. 2001; Turner et al. 1999a; Ballantyne et al. 2001). However, the ionized Fe Ka line 
is not unique to NLS1 galaxies, and some BLS1 galaxies also show ionized Fe Ka. This 
may suggest that the luminosity of the central source plays as important a role as the 
accretion rate (e.g., Guainazzi et al. 1998). 
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Alternative explanations for the extreme properties of NLSls are: (i) the size of the 
BLR of NLSls is larger (i.e., the BLR gas is more distant from the nucleus) than that in 
the BLSls (Guilbcrt, Fabian, & McCray 1983; Mason et al. 1996; Wandcl & Boiler 1998) 
resulting in the narrowness of the width of the permitted lines due to a lower orbital 
velocity; (ii) we have a nearly face-on view of a flattened BLR in NLSls (Osterbrock 
& Pogge 1985). Assuming the motion of the BLR gas around the central super-massive 
black-hole to be virialized, the narrowness of the lines is due to the fact that the gas is 
moving preferentially on a plane that is almost perpendicular to the line of sight resulting 
in the smaller velocity dispersions, hence the line widths are reduced by a factor sim, 
where i = is face-on. However, Boroson & Green (1992) and Kuraszkicwcz et al. (2000) 
do not favor the low inclination model, while Nandra et al. (1997) showed that the 
inner regions of BLSls also appear to be observed nearly face-on. Reverberation results 
(Kaspi et al. 2000; Peterson et al. 2000) find that the BLRs of NLSls and BLSls have 
comparable sizes, while NLSls have virial masses that are one order of magnitude smaller 
than BLSls. This result shows that the size of the BLR does not scale with the mass 
of the central black-hole but with the luminosity (Laor 1998) which is connected to the 
accretion rate. Dewangan et al. (2001a, 2001b) suggested that both the steeper X-ray 
emission and narrower width of the H/3 line, and also the anti-correlation between the 
slope of the X-ray spectrum and the width of the H/3 line could be due to the variation 
in the fractional accretion rate. In this scenario, the higher fractional accretion rate 
producing an increased accretion disk emission in the soft X-rays and a steeper power 
law in the hard X-rays, results in a higher radiative pressure per unit gravitating mass. 
Thus, for a given density distribution of diffuse material, the BLR clouds will form further 
out in a higher-intensity radiation field, resulting in the narrower emission lines. 

Variability studies using long ASCA (Advanced Satellite for Cosmology and 
Astrophysics) observations have been a valuable tool to explore the nature of the soft 
X-ray excess emission and other spectral components of NLSls (Turner et al. 2001b; 
Romano et al. 2002). A 35-day long ASCA observation of a NLS1 galaxy, Akn 564, re- 
vealed that a slower varying soft excess component is superimposed on a fast varying 
continuum component (Turner et al. 2001b). Similar results are also inferred from a 
12-day ASCA observation of another NLS1 galaxy Ton S180 (Romano et al. 2002). 

IRAS 13224-3809 is an extremely variable NLS1 galaxy at a redshift of 0.06667 and 
with soft (0.1 - 2.4 keV) X-ray luminosity of 3 x 10 44 erg s' 1 (Boiler et al. 1993). The 
FWHM of the H/3 line of this source is only ~ 700 km s _1 (Boiler et al. 1993; Leighly 
1999b) which is comparable to the width of the forbidden line [O III]A5007. ROSAT ob- 
servations (Boiler et al. 1993) showed a complex soft X-ray emission that was very steep 
(rx ~ 4.4) and rapidly variable (change in intensity by a factor of 2 in ~ 800 s ). 
Subsequent ASCA observations in 1994 confirmed the complex soft X-ray emission and 
variability (Leighly et al. 1997; Leighly 1999a,1999b). A 30-day ROSAT RRl monitoring 
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of IRAS 13224-3809 has revealed the most extreme and multiple giant amplitude X-ray 
variability (Boiler et al. 1997). 

In this paper we present the results from a 10-day ASCA observation of the NLS1 
galaxy IRAS 13224-3809. In Sect. || we describe our observations and data reduction. In 
Sect. [| we discuss the time variability of the source. We analyze the mean spectrum in 
Sect. ^ and time resolved spectra in Sect. || In Sect. ^ we present flux resolved spec- 
troscopy. In Sect. [?] we compare our results for IRAS 13224-3809 with the properties of 
Akn 564 and Ton SI 80. Finally, we discuss the results in Sect. || and summarize our 
results in Sect. [| 

2. Observation & Data Reduction 

ASCA consists of four focal plane detectors, two CCDs (the solid-state Imaging 
Spectrometers, SISO and SIS1, 0.4 - 10 keV, Bruke et al. 1991) and two GISs (the Gas 
Imaging Spectrometers, GIS2 and GIS3, 0.7— 10 keV, Ohashi et al. 1996, and references 
therein). All the four detectors operate simultaneously. ASCA observed IRAS 13224-3809 
(Principal Investigator: K. M. Leighly) for a total duration of ~ 834 ks starting from 
JD=2451731.562 (for the screened data) in the 1CCD mode. The data were reduced us- 
ing standard techniques (Revision 2). Data screening yielded an effective exposure time 
of 228 ks for SISO, 226 ks for SIS1, and 271 ks for both the GISs. The mean SISO count 
rate was (5.063±0.061) x 10 -2 count s _1 , which is about 30% higher than the SISO count 
rate of (3.051 ± 0.095) x 10~ 2 count s^ 1 found during a previous ASCA observation in 
1994 (Leighly et. al. 1997; Leighly 1999a, 1999b). 

After the launch in 1993, ASCA SIS detectors degraded gradually in efficiency at 
lower energies, due to the increased dark current levels and charge transfer inefficiency 
(CTI). This degradation resulted in SIS spectra which diverge from each other and from 
the GIS data. The instruments can diverge by as much as 40% for energies below 0.6 keV 
for data taken in 2000 Januarj^]. The degradation in efficiency is not well understood 
and it could not be corrected for by any of the software at the time of writing this 
paper. There has been a non-linear evolution of the SIS CTI during the last phase of 
ASCA observations (AO-8). The SIS team has revised the calibration of the non- uniform 
CTI effect and released a new calibration file (sisph2pi_l30201.fits) on 2001 March 29. 
The IRAS 13224-3809 data were calibrated using the above revised calibration file. 

The divergence of the SIS detectors at low energies can be compensated for in the 
spectral analysis by employing the technique of Yaqoob et al. (2000), who provide an 
empirical correction by parameterizing the efficiency loss as a time-dependent absorption 
term0 ("excess Nh")- The correction for SISO follows a linear relationship, Ah (SISO) = 



1 see 


http://heasarc.gsfc.iiasa.gov/docs/asca/watchout.htmj 


2 see 


http://lheawww.gsfc.nasa.gov/~yaqoob/ccd/nhparam.html 
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Fig. 1. SIS0+SIS1 light curve (in 500 s bins) of IRAS 13224-3809 obtained from the 
ASCA data taken during 2000 July 6-15. The light curves were extracted in the energy 
band of 0.7 — 10 keV. The horizontal axis shows the time in seconds from the beginning 
of the observation (MJD = 51731.188273 day). 

3.635857508 x 10" 8 ( T - 3.0 1 74 8 28 x 10 7 ) 10 20 cm~ 2 , where T is the average of start 
and stop times of the observation measured in seconds since launch. The SIS1 excess 
absorption term does not follow the simple linear relationship with time but it is found 
that a slightly larger absorption column can be applied to the SIS1 data so that both 
the SIS detectors agree well at lower energies. For the observations of IRAS 13224-3809, 
TVh(SISO) = 7.53 x 10 20 cm" 2 , where T = 2.37 x 10 8 s and we adopted 7V H (SIS1) = 
1.05 x 10 21 cm- 2 . 

3. Time Variability 

Light curves were extracted using bin sizes of 500 s in the band (0.7 — 10 keV) from 
the SIS data and 5000 s in the soft (0.7 - 1.3 keV) and hard (1.3 - 10 keV) bands from 
both the SIS and GIS data. The soft and hard bands were chosen to have similar signal- 
to-noise and to separate approximately the two spectral components - soft excess and 
power law (see Sect. ^|). The exposure requirements for the light curves were that the 
bins be at least 50% and 10% exposed in each instrument for the 500 s and 5000 s 
curves, respectively. Background light curves were extracted from the source free regions 
and subtracted from the source light curves after appropriate scaling to compensate for 
different sizes of extraction regions. We combined the 500 s light curves from the two SIS 
detectors only, and 5000 s curves in each band from all the four detectors. The observed 
counts correspond to a mean observed flux of 6.7 x 10 -13 erg s _1 cm -2 , and luminosity 
of 5.9 x 10 42 erg s^ 1 (assuming H Q = 75 km s _1 Mpc -1 , q = 0.5) in the 2-10 keV 
band. 

Fig. [I] shows the background subtracted 0.7 — 10 keV band light curve with 500 s 
bins. Fig. || shows the background subtracted SIS soft-band and SIS+GIS hard-band 
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Fig. 2. Light curves of IRAS 13224-3809 with 5000 s bins. The top panel is the SIS 
soft-band (0.7 — 1.3 keV) light curve; the second panel is the SIS+GIS hard-band (1.3 — 
10 keV) light curve; the third panel is the total band (0.7 — 10 keV) light curve. The 
last panel shows the hardness ratio defined as the ratio of the count rates in the 1.3 — 
10 keV and 0.7— 1.3 keV bands. The light curves have been corrected for the background 
contribution. The horizontal axis shows the time in seconds from the beginning of the 
ASCA observation (MJD = 51731.188273). 



light curves in 5000 s bins. Also shown in Fig. || is the 0.7 — 10 keV band light curve, 
and hardness ratio (HR) defined as the ratio of count rates in the 1.3 — 10 keV and 
0.7- 1.3 keV bands. 

The light curves binned to 5000 s (Fig. ||) show trough-to-peak variations in the 
count rate by a factor of 25 in the soft band, ~ 20 in the hard band. The light curve 
sampled on 500 s (Fig. |l|) reveals even higher amplitude variations due to fast flickering, 
with a maximum amplitude of variability of a factor J> 37. Close up examination of the 
light curve in Fig. |l| reveals several rapid variable events. The most rapid event, shown in 
Fig. ||, occurred at 2.755 x 10 5 s after the beginning of the observation. The count rate 
increased from 0.12 to 0.24 count s _1 just within 500 s . Several variable events with a 
change in the count rate by a factor of ~ 3 on a timescale of 3000 — 5000 s have been 
detected e.g., events seen at 3.265 x 10 5 s , 3.69 x 10 5 s , and 5.775 x 10 5 s . The soft 
and hard X-ray light curves, shown in Fig. [2] , show similar variability properties. We 
have calculated the cross correlation function (CCF) of the hard X-ray (1.3 — 10 keV) 
flux with respect to the soft X-ray (0.7 - 1.3 keV) flux. The CCF is plotted in Fig. |, 
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Fig. 3. The most rapid variability event observed during the 10-day ASCA observation 
of IRAS 13224-3809. The light curve is an expanded view of the light curve shown in 
Fig. [l] during the variability event. 
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Fig. 4. Cross-correlation function of hard (1.3 — 10 keV) flux with respect to the soft 
(0.7 — 1.3 keV) flux observed during the year 2000. 



which shows strong correlation between the hard and soft X-ray fluxes without any 
significant time delay. Since the power-law component contributes ~ 32% of the total 
flux in the soft (0.7 — 1.3 keV) band as inferred from the mean spectrum (see Sect. 
the observed correlation is partly due to variation in the power-law component alone. 
However, variability of the power-law flux alone is not sufficient to explain the soft X-ray 
variability. This can be seen from Fig. which shows an expanded view of the flaring 
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Fig. 5. An X-ray flare, in the soft and hard bands, from IRAS 13224-3809 observed 
1.9 x 10 5 s after the beginning of the observation in the year 2000. 
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Fig. 6. Same as Fig. except that the light curves are shown with time bins of 1 day. 
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event seen at 1.9 x 10 5 s after the beginning of the observation. The 0.7 — 1.3 keV band 
flux changed by a factor of ~ 4 from 0.010 ± 0.003 count s -1 to 0.044 ± 0.006 count s -1 , 
while the 1.3 — 10 keV band flux changed by a factor of 6 from 0.016 ±0.005 count s -1 to 
0.097 ± 0.009 count s _1 . If the soft X-ray flux above the hard X-ray power law remains 
constant, so that the observed soft X-ray variability is entirely due to changes in the 
power-law flux alone, then a factor of ~ 12 change is required in the power-law flux in 
the 0.7—1.3 keV band. This required change is much higher than the factor of 6 observed 
in the 1.3 — 10 keV flux. Therefore, the soft-excess flux above the hard X-ray power law 
must also have changed by a factor of ~ 2.6 either simultaneously or with a short time 
delay with respect to the hard 0.3 — 10 keV flux. 

From Fig. || and Fig. [|, it is clear that some events have a sharper rise in hard X-rays 
e.g., the flaring event seen at 1.9 x 10 5 s after the beginning of the observations, with 
accompanying change in the hardness ratio, while during other flaring events between 
3 — 6 x 10 5 s , there is no change in the hardness ratio. When the light curves are binned 
on a timescale of a day, another important type of behavior of the source is observed. 
Fig. ^| shows the light curves and HR with time bins of a day. The HR appears to increase 
when the flux is rising or falling but settles down to a lower value when the maximum 
flux is reached. The hardness ratio tells only the relative changes in the soft and hard 
bands. Any possible change in the spectral shape will be investigated in Sect. 

3.1. Fractional Variability Amplitude 

The fractional variability amplitude F vav and its error OF var are defined as 



(Edelson et al. 2001) where S 2 is the total variance of the light curve, (a 2 lr ) is the mean 
squared error, and {X) is the mean count rate. 

First we calculated the fractional variability amplitude of the total band (0.7— 10 keV) 
light curve with 500 s bins to be F var = 73.5 ± 2.6%. This quantity measures deviations 
relative to the mean, integrated over the entire duration of the observation. We also 
calculated F val in the soft (0.7 — 1.3 keV) and the hard (1.3 — 10 keV) bands, from the 
light curves shown in Fig. || with 5000 s bins. F var thus calculated is 75.3 ± 4.6% in the 
0.7 - 1.3 keV band, and 66.1 ± 3.8% in the 1.3 - 10 keV band. 

We also measured F var for each day after splitting the light curves with 5000 s bins 
into 10 evenly-sampled sections across the 10-day ASCA observation. Fig. ^ shows F var 
in the soft, and hard bands calculated for each day. The soft-band F var changes across 
the 10-day observation, a constant fit resulting in a minimum x 2 of 21.88 for 9 dof. The 
hard-band F var , however, does not appear to vary significantly, a constant fit to the hard- 
band .Fvar curve giving a minimum \ 2 of 9-13 for 9 dof. However, there appears to be 
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Fig. 7. Variability properties of IRAS 13224-3809. First two panels from the top show 
fractional variability amplitude in different energy bands as a function of time. The last 
panel shows the SISO count rate in the 0.7 — 10 keV band. 

some similarity in the soft and hard-band F val curves in Fig. ^. As already mentioned, 
the power-law component contributes significantly to the unabsorbed soft-band X-ray 
flux (see Sect. ^). Therefore, the gross similarity between the soft-band and hard-band 
F var and a stronger variability of the soft-band F var can be understood if the changes 
in the soft-band F var are not only caused by variations in the power-law continuum flux 
but also by the intrinsic variations in the soft-excess component above the power law. It 
is quite possible that the variations in the above two spectral components are correlated 
probably with a short time delay. The quantity F var in the soft or hard band, is not 
correlated with the observed count rate as can be seen in Fig. suggesting that the 
variability properties do not depend on the flux level. 

4. The Mean Spectrum 

For each detector, photon energy spectrum of IRAS 13224-3809 was accumulated from the 
entire observation. Pulse invariant channels were appropriately grouped for the spectral 
analysis while considering the degradation in the energy resolution of the SIS detectors. 
The data from the four instruments were fit simultaneously while keeping the relative 
normalizations free allowing for the small differences in the calibration of the absolute 
flux, and differences in the fraction of encircled counts contained in the SIS and GIS 
extraction cells. The spectral fits were performed with the XSPEC VI 1 .0. 1 package, us- 
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Fig. 8. Deviation x of the mean spectrum of 2000 July from the best-fit 2 — 10 keV power 
law modified by Galactic absorption where \ = — i — — i — j N° is the observed counts 
in energy channel i, Ci is the standard error on A^ obs , and _/V" lod is the best-fit model 
counts in channel i. 



ing response matrices released in 1997 for the GISs, and response files generated using 
HEAsoft v5. 0.4 for the SISs. 

The spectral shape was first determined by fitting a redshifted power-law model mod- 
ified by Galactic absorption (N H = 4.79 x 10 20 cm -2 ; Dickey & Lockman 1990; Model 
A) to the data above 2 keV. An additional absorption term was used for the SIS0 
and SIS1 to compensate for the low energy degradation as described in Sect. ||. The 
models for Galactic absorption use the absorption cross-sections of Balucinska-Church 
& McCammon (1992). For this excercise we used the SIS data in the energy band of 
2 — 7.32 keV and GIS data in the 2 — 10 keV band (both in the observer's frame). The 
power-law fit yielded T x = 2.02 ± 0.08 and a minimum \ 2 of 243.52 for 183 dof. The 
errors quoted, here and below, were calculated for the 90% confidence level based on 
y 2 • +2.71. 

The results of the fit described above are given in Table g for Model A. We list the 
best-fitting power-law photon index (rx) in the 2 — 10 keV band, the unabsorbed power- 
law flux (/x L ) m t ne ener gy band of 2 — 10 keV, the minimum reduced x 2 (Xred)> ano - °-°f- 
The deviation of the best-fit model from the observed data is shown in Fig. || in terms of 
X = — ' — ^-^ — , where N° is the observed counts in energy channel i, <ii is the standard 
error on N° hs , and jV z mod is the best- fit model counts in channel i. The data below 2 keV 
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Table 1. Best-fit model spectral parameters for IRAS 13224-3809 derived from the 10- 
day ASCA observation. 
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1 Model A is the best-fit redshifted power-law model in the 2 — 10 keV band modified by 
the Galactic absorption. Model B is the combination of redshifted blackbody and power law 
modified by the Galactic absorption. The best-fit parameters for the model B were derived 
using the total energy band of 0.7 — 10 keV. 



2 Intrinsic flux in the energy band of 0.7 — 10 keV and in the units of 10 13 erg cm 2 s 1 . 

3 Intrinsic luminosity in the rest frame and in the energy band of 0.7 — 10 keV and in the units 
of 10 42 erg s" 1 . 

4 Intrinsic flux of the soft hump (described by a blackbody) in the 0.7 — 2 keV band and in the 
units of 10~ 13 erg cm~ 2 s _1 . 

5 Intrinsic flux of the power-law component in the 2—10 keV band and in the units of 
m-13 „„„ ™-2 -i 
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Fig. 9. Mean spectra of IRAS 13224-3809 obtained in the year 2000 and the best-fit 
model - blackbody and power law modified by the Galactic absorption (Model B) fitted 
over the entire energy range (top panel) and deviation of the observed mean spectrum 
from the best-fit model (bottom panel). 

is also shown in Fig. |^. This plot helps to indicate the significant features in the spectrum. 
A strong soft excess is evident, appearing as a hump of emission below 2 keV. This feature 
was also found in the ROSATPSPC observations (Boiler et al. 1993), and in the previous 
ASCA observation of 1994 (sequence number 2011000, Leighly 1999b). Hereafter, we refer 
to this component as the "soft hump" . Minor calibration problems are also visible mainly 
below 2 keV. We do not detect an Fe Ka line from IRAS 13224-3809 in the energy range 
6.4 — 7 keV where the deviations are below 3cr level (see Fig. ||). The 90% confidence 
upper limit for the equivalent width of Fe Ka line is found to be ~ 400 eV. 

4.1. The Soft Hump 

We confirm the presence of the soft X-ray excess emission component previously observed 
by Boiler et al. (1993) and Leighly (1999b). Similar features are well known in other NLS1 
galaxies, e.g. in RE J1034+393 (Pounds et al. 1995), Akn 564 (Turner et al. 2001b), and 
Ton S180 (Romano et al. 2002). In these objects, Chandra LETG results (Ton S180, 
Turner et al. 2001a; NGC 4051, Collinge et al. 2001) show that the soft hump is a 
smooth continuum component, as opposed to a blend of unresolved spectral features. 
This allows us to choose a continuum component model to parameterize the soft hump 
emission. We use the blackbody model to characterize the soft hump component as it 
adequately models the shape and flux of the soft hump. We used the SIS data in the 
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range 0.7 — 7.32 keV simultaneously with the GIS data in the range 1 — 10 keV, and 
fitted the redshifted blackbody and power-law model modified by the Galactic absorption 
(Model B). An additional absorption term as described above was also used for the 
SIS data. The best-fit blackbody and power-law model yielded a rest-frame temperature 
kT = 130.713 2 e V> absorption corrected blackbody flux, f x — 4.34 xl0~ 13 erg s _1 cm~ 2 
in the energy band of 0.7-1.3 keV, and T x = 2.11±0.05 for minimum \ 2 = 421.98 for 244 
dof. The observed data and the best-fit model are shown in Fig. |^. Also shown in Fig. || 
are the deviations of the observed data from the best-fit model. The fit is poor mostly due 
to calibration problems at low energies and residuals near 1 keV. The absorption feature 
near 1 keV has been detected in the earlier ASCA observation of 1994 and has been 
interpreted as the blueshifted absorption edges of oxygen (see Leighly et al. 1997). Here 
we note that the GIS and SIS detectors do not agree at the position of the absorption 
feature and we do not fit an absorption line or edge model. The uncertainty in the low 
energy calibration and the degradation in the energy resolution of the SIS detectors make 
it difficult to explore the absorption feature. We also note that the parameterization of 
the soft hump as a blackbody does not alter the power-law slope significantly. We find 
that the power-law continuum contributes 32.3% of the flux in the soft 0.7 — 1.3 keV 
band in the mean spectrum while the soft hump above the power-law contributes only 
2.7% in the 1.3 — 10 keV band. Although the power law has a significant contribution to 
the soft X-ray flux, the blackbody contribution to the the hard X-ray flux is negligible. 

5. Time Resolved Spectroscopy 

5.1. Method and Selection Details 

To examine the spectral variations of IRAS 13224-3809, we extracted 10 time-selected 
spectra across the 9.6-day ASCA observation using Xselect 2.0. As the source is not 
bright enough, it was not possible to extract spectra following individual flares and dips 
with sufficient signal-to-noise ratio. Instead, we chose the sampling timescale of one day, 
except for the last day of observation for which the sampling time was 0.6 day. The 
vertical dotted lines in Fig. [l] show our 10 intervals within which spectra were extracted. 
The resulting "on-source" average exposure time was ~ 20 ks per detector. We set the 
ancillary response, and response matrix files to be those of the mean spectrum. Again, 
the energy channels were appropriately grouped to achieve a good signal-to-noise (a 
minimum of 20 counts per energy channel) while considering the spectral resolution of 
each detector. All fits were performed by fixing the relative instrument normalizations 
obtained from the best-fit values from the mean spectrum. The correction for the low 
energy SIS problem was performed in the same way as for the mean spectrum. All models 
included the Galactic absorption as before. Fig. [l(] shows the results of our analysis in 
the form of a time series for the various parameters and are described in detail below. 
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Fig. 10. Spectral variability of IRAS 13224-3809. From the top, the time series are the 
photon index, power-law flux in units of 10 -13 erg cm~ 2 s -1 and in the 2 — 10 keV band, 
blackbody temperature in eV, blackbody flux in units of 10~ 13 erg cm~ 2 s _1 and in the 
0.7 — 2 keV band, and the total flux in the 0.7 — 10 keV band. All the fluxes have been 
corrected for the Galactic absorption. 

Time assignments of the spectra refer to the mid point of the observation, in days, from 
the beginning of the observation. 

5.2. Variability of the Continuum 

We fitted each of the 10 time-selected spectra by a redshifted power-law model (Model 
A) using the 2-7.32 keV band for the SIS data and the 2-10 keV band for the GIS data 
(the same data exclusion as for the mean fit). The results of the fits are listed in Table [j] 
. In order to explore the time evolution of the soft hump, we also fitted each spectrum 
by a combination of redshifted blackbody and power-law models (Model B). For these 
fits, SIS data in the band 0.7 — 7.32 keV and GIS data in the band 1 — 10 keV were 
used for each time-selected spectrum. The results of these fit are also listed in Table [l] . 
As can be seen in Table [l], there are no significant changes in the best-fit photon indices 
obtained from fitting models A and B. Fig. |lO| shows the time series for the photon index, 
and model power-law flux in the 2 — 10 keV band, blackbody temperature, blackbody 
flux in the 0.7 - 2 keV band, and total flux in the 0.7 - 10 keV band. Note that the 
best-fit values for the power-law model plotted in Fig. [l^ are those derived from model 
B. The best-fit values of Fx range from 1.74 to 2.47 across the 10-day observation, but 
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Table 2. Matrix of linear correlation coefficients and the corresponding significance level 
calculated for the best-fit parameters obtained from the 10 time-selected spectra. 
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the day-to-day variations are not significant. A significant variation in the Tx is observed 
between day 5 and day 9, the change in Fx being 0.57 ± 0.21. We also note that the 
power-law component dominates the 2 — 10 keV band and the flux variations in this 
band on a timescale of ~ a day are due to the changes in the continuum level. 

5.3. Variability of the Soft X-ray Hump 

The blackbody flux (/x B ) m the 0.7—2 keV band roughly follows the power-law flux (/x L ) 
in the 2 — 10 keV band (see Fig. [b]). The trough-to-peak variation in the soft hump flux 
is by a factor of ~ 4.7, while that of the power-law flux is by a factor of ~ 3.5 suggesting 
that the soft hump is more variable than the power-law component on timescales of 
~ a week. The blackbody temperature varies from 98.4 eV to 137.6 eV across the 10- 
day observation, but these variations are not significant considering the error bars. The 
contribution of the power-law component to the total flux in the 0.7 — 1.3 keV band 
varies from 24% on day 10 to 48% on day 3 and is not correlated with any of the other 
spectral parameters. 

To examine the time evolution of soft hump flux, we constructed a plot to highlight 
the variation of the soft hump above the power law. Fig. |ll| shows the deviations of the 
observed data from the best-fit 2 — 10 keV power-law model for all the 10 spectra. The 
data below 2 keV have also been plotted. We note that the soft hump is always evident 
above the power-law continuum. Strong variations in the soft hump are also evident. 

5.4. Correlation between spectral parameters 

In order to investigate possible correlation between spectral parameters, we have cal- 
culated linear correlation coefficients. Table ^ shows the matrix of linear correlation 
coefficients and the corresponding significance levels. Fig. |l2] shows plots of spectral pa- 
rameters. The most significant correlations are those involving fluxes, both soft hump 
flux and power-law flux are strongly correlated with the total flux. The power-law pho- 
ton index is better correlated with the soft hump flux (at 96.6% level) than with the 
power-law flux (at 91.8% level) suggesting that the shape of the power law is probably 
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Fig. 11. Deviations \ of the time-selected spectra from the 2 — 10 keV best- fit power-law 
model modified by the Galactic absorption. A strong variation in the soft hump emission 
above the power law is evident. 

more sensitive to the seed photons than to the changes in the power-law flux. The soft 
hump flux is correlated with the blackbody temperature as is expected for blackbody 
emission. 
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Fig. 12. Correlations between spectral parameters: (a) 2 — 10 keV photon index (Fx) 
vs 2 - 10 keV power-law flux (/| L ), (b) T x vs blackbody flux in the 0.7 - 2 keV band 
(/x B )i ( c ) r x vs total flux (blackbody + power law) in the 0.7 - 10 keV band (/ x otal ), 
(d) /P L vs /| B , (e) / x otal vs /| B , and (f) / x otal vs /^ L . 

6. Flux Resolved Spectra 

The soft hump flux and the power-law flux appear to be correlated with the photon 
index at significance levels of 96.6% and 91.8% (Table 2; see also Fig. |l2|). In order to 
improve the statistical significance and to further explore the dependence of spectral 
parameters, we have carried out spectral analysis at different flux levels. We extracted 
three averaged spectra corresponding to the SIS0 count rates of < 0.07 count s _1 (low), 
0.07 — 0.14 count s" 1 (intermediate), and > 0.14 count s _1 (high). The spectra were 
analyzed in the same way as before and the results of the spectral fitting are given in 
Table [l]. Fig. [l3] shows the deviation of the data from the best-fit 2 — 10 keV power law. 
In all the three states, the soft hump is present above the power law and is strongly 
variable. There is a significant change in the power-law slope between the low and the 
high states, Ar x = 0.44 ± 0.12. The soft hump flux in the 0.7 - 2 keV band also varied 
between the two states by a factor of ~ 5.6, while the power-law flux in the 2 — 10 keV 
band varied by a factor of ~ 3.2. Fig. |l3| also shows the contours of allowed values of Tx 
and blackbody normalization at 68%, 90%, and 99% confidence levels. It is clear from 
the contour plot that when the intensity of the source increases, the power law becomes 
steeper and the soft hump stronger. 
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Fig. 13. Deviations x of the flux-selected spectra from the best-fit 2—10 keV power law 
in the low, intermediate, and high state and the confidence contours at 68%, 90%, and 
98% levels of the photon indices and blackbody normalizations derived from the best-fit 
Model B. 

Table 3. Best-fit model spectral parameters for IRAS 13224-3809 derived from the 1994 
ASCA observation. 
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1 Model A is the best-fit redshifted power-law model in the 2 — 10 keV band modified by 
the Galactic absorption. Model B is the combination of redshifted blackbody and power law 
modified by the Galactic absorption. The best-fit parameters for the model B were derived 
using the total energy band of 0.7 — 10 keV. 



2 Intrinsic flux in the energy band of 0.7 — 10 keV and in the units of 10 13 erg cm 2 s 1 . 

3 Intrinsic luminosity in the rest frame and in the energy band of 0.7 — 10 keV and in the units 
of 10 42 erg s _1 . 

4 Intrinsic flux of the soft hump (described by a blackbody) in the 0.7 — 2 keV band and in the 
units of 10~ 13 erg cm~ 2 s _1 . 

5 Intrinsic flux of the power-law component in the 2 — 10 keV band and in the units of 
10~ 13 erg cm~ 2 s~\ 

6 Reduced minimum \ 2 /degrees of freedom. 
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Fig. 14. Light curve of IRAS 13224-3809 sampled with 500 s and in the total energy 
bands derived from the observation of 1994. The vertical dotted lines show the time 
selection for the low and high flux states. 



In order to further confirm the above results, we have also analyzed the data obtained 
from an earlier observation in 1994 with ASCA . First we extracted light curves using 
bin sizes of 500 s in the total band (0.7 — 10 keV) for the SISs. Background light curves 
were also extracted and the source light curves have been corrected for the background 
contribution. The exposure requirements were that the 500 s bins be at least 50% exposed. 
The final light curve was constructed after combining the light curves from the two SIS 
detectors and is shown in Fig. The timing properties of the source using these data 
have been studied in detail by Leighly (1999a). Our aim here is to show our time selection 
for the low and the high intensity states of the source used for spectral analysis. We 
extracted spectra from the two time intervals where the count rate is low (low state) 
and also from the time interval where the source count rate is high (high state). We 
also extracted spectra from the total time interval of the observation (mean spectrum). 
The spectra were analyzed in the same way as before except for the SIS gain correction 
which was not used. The mean, low and high state spectra were fitted by a redshifted 
power law modified by the Galactic absorption and in the energy band of 2-10 keV(Model 
A). The best- fit parameters are listed in Table ||. The deviations of the observed data 
from the best-fit power law are shown in Fig. [l5|. The soft excess, seen in Fig. [ll], was 
again parameterized by a redshifted blackbody. The best-fit parameters obtained from 
the blackbody and power law modified by the Galactic observation (model B) are also 
listed in Table |[ It is evident from Table || that the photon index becomes steeper with 
corresponding increase in the blackbody flux in the high state compared to the low state. 
The same result is also seen in the \ 2 contour plots of Tx and blackbody normalization 
in Fig. [l5|. The above result is similar to the results shown in Fig. [l3| 
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Fig. 15. Deviation \ of the spectra, observed in 1994, from the best-fit 2 — 10 keV power 
law. See text and Fig. [l4| for definition of different spectra. Also shown are the confidence 
contours at 68%, 90%, and 98% levels of the photon index and blackbody normalization 
from the blackbody and power-law model fits to the different spectra. 

7. Comparison with Akn 564 and Ton S180 

A comparison of our results for IRAS 13224-3809 with those obtained for Akn 564 (Turner 
et al. 2001b) and Ton S180 (Romano et al. 2002) reveals a broad similarity in the overall 
shape of the X-ray spectra (strong soft X-ray excess, steep power law) and variability 
properties. However, there are some important differences: (i) An Fe Ka line is not 
detected from IRAS 13224-3809, while both Akn 564 and Ton S180 show the presence 
of an Fe Ka line with a large equivalent width from highly ionized material, (ii) The 
contribution of the power law to the soft hump emission is only 32% in IRAS 13224- 
3809 while it is 75% for Akn 564 and 72% for Ton S180. (Hi) The mean photon index 
(Ix = 2.11±g;g|) of IRAS 13224-3809 is flatter than that of Akn 564 (T x = 2.538±0.005) 
and Ton S180 (Fx = 2.44 ±0.02). (iv) The power-law slope appears to be correlated with 
the soft hump flux in Ton S180 and IRAS 13224-3809 while the correlation is absent in 
Akn 564. (v) The variability amplitude of the soft hump and power-law components are 
higher in IRAS 13224-3809 than that of Akn 564 and Ton S180. 

8. Discussion 

IRAS 13224-3809 shows the rapid and large amplitude variability seen earlier with 
ROSAT (Boiler et al. 1993; Boiler et al. 1997) and ASCA (Leighly 1999a). The 
0.7 — 10 keV band light curve with 500 s bins shows trough-to-peak variation by a 
factor 37 during our 10-day ASCA observation. Rapid X-ray variability by a factor 
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of 2 — 3 on a timescalc of ~ 2000 s has also been observed. During the 10-day obser- 
vation, the light curves sampled with 5000 s bins show trough-to-peak variations by a 
factor i£ 25 in the soft band (0.7 — 1.3 kcV), and about a factor of 20 in the hard band 
(1.3 — 10 keV). Variability events appear to be sharper in the hard X-ray band than in 
the soft X-ray band, and the intensities in the two bands are strongly correlated. Changes 
in the power-law component alone are not sufficient to produce the observed correlation 
but the soft hump above the power law changes its flux simultaneously or with a short 
time delay (see Fig. || and Sect. ||). 

The mean photon index of the X-ray power law, obtained from the spectral fit to 
the mean spectrum, is Fx = 2.1lj^Qg. The contribution of the hard X-ray power-law 
component to the flux in the soft band (0.7 — 1.3 keV) is only about 32% which is much 
smaller than that found for Akn 564 (75%, Turner et al. 2001b) and Ton S180 (72%, 
Romano et al. 2002), suggesting that the soft hump component is more pronounced in 
IRAS 13224-3809. Our time resolved spectroscopy reveals variations in Tx from 1.74+q'^| 
to 2.47^q'^3, implying a Arx = 0.73 ± 0.37. The mean power law appears to be flatter 
than that obtained for other NLS1 galaxies like Akn 564 (T x = 2.538±0.005, Ar x = 0.27, 
Turner et al. 2001b) and Ton S180 (T x = 2.44 ± 0.02, Ar x = 0.24, Romano et al. 2002) 
and the variation in Tx is slightly higher in spite of the fact that the time span of the 
ASCA observations for the later two objects were longer. Thus IRAS 13224-3809 shows 
higher amplitude variability in the power-law slope as well as in the soft hump and the 
power-law intensity. This indicates that the physical parameters governing the X-ray 
emission vary by larger factors in IRAS 13224-3809 compared to that for Akn 564 and 
Ton S180. 

8.1. The Fx - Luminosity Relation 

Our time resolved spectroscopy has revealed that the 2 — 10 keV continuum steepens 
with increase in the flux. The photon index changes by ~ 0.4 with the corresponding 
change in the power-law flux being a factor of ~ 3.2 between the low and high flux states 
observed in the year 2000. Similar behavior is seen during the 1994 observations. Fig. |l6| 
shows the plot of Tx against the 2 — 10 keV flux. The steepening of Tx with flux has 
also been observed in a number of Seyfert 1 galaxies (e.g., Singh et al. 1991; Done et 
al. 2000; Zdziarski & Grandi 2001; Petrucci et al. 2001; Vaughan & Edelson 2001, and 
Nandra 2001) and can be understood in the framework of thundercloud and accretion 
disk model of Merloni & Fabian (2001). The basic building blocks of this model are the 
active regions above an accretion disk, viewed as magnetic thunderclouds and triggered 
by magnetic reconnection. The sizes of the active regions are distributed as a power law. 
Rapid X-ray flares are produced in the active regions by inverse Compton scattering of 
soft photons from the disk and Tx-luminosity relation of the form Tx = To — KL~ S 
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is expected (Merloni & Fabian 2001), where the asymptotic value of Tx, To, depends 
mainly on the optical depth tt of the active regions and disk or seed photon intensity. 
The exponent S, which determines the amount of spectral variation, is mainly dependent 
on the spatial distribution of correlated flares. We have fitted the above Tx~L relation to 
the IRAS 13224-3809 data (see Fig. ^ and obtained the best-fit values T = 3.24 ±0.97, 
6 = 0.39 ± 0.32, where the quoted errors are at la level. Thus for IRAS 13224-3809, 
the asymptotic photon index is similar within errors to that obtained for the Seyfert 1 
galaxy MCG-6-30-15 (r = 2.30±g;^; Merloni & Fabian 2001). Given the large error 
bars in T for IRAS 13224-3809 as well as for MCG-6-30-15, the coronal optical depth 
for IRAS 13224-3809 does not appear to be significantly different from that inferred 
for MCG-6-30-15 {tt > 1-5). However, better quality data such as that obtained from 
monitoring observations with XMM-Newton are required to make a firm conclusion. In 
the framework of thundercloud model a smaller covering fraction of the active regions is 
required for IRAS 13224-3809 than for MCG-6-30-15. The smaller the covering fraction, 
the larger the observed variability and greater the chance of a large flare to occur (Merloni 
& Fabian 2001). 

8.2. The Soft Hump Emission 

Our time resolved spectroscopy has revealed that the soft hump emission and the power- 
law flux are correlated (see Figs. [nj [l3], [l^). Pounds et al. (1995) first noted that the 
reprocessing of hard X-rays is insufficient to produce the observed soft hump emission of 
NLS1 galaxies and the soft hump emission could be the intrinsic disk emission resulting 
from near or super-Eddington accretion rates. In this scenario, sharper variability events 
are expected in the soft X-rays than in the hard X-rays which is contrary to that observed 
from IRAS 13224-3809 (see Fig. ||). On the other hand, in the framework of reprocessing 
models, sharper variability events are expected in the hard X-rays than in the soft X- 
rays. Although a correlation between the soft and hard X-ray flux is expected in both 
the above scenarios, it is difficult to understand the origin of soft X-ray emission either 
as the intrinsic disk emission or as the reprocessing of hard X-rays. 

We propose yet another mechanism which may be partly responsible for the observed 
soft hump emission from NLS1 galaxies. If the accretion rate is super-Eddington, the 
accretion flow is likely to be dense and optically thick. The inner regions of such a disk 
are supported by the strong radiation pressure. Consequently, the disk puffs up and 
becomes geometrically thick (see e.g. Collin 2001). The electron density of the accretion 
disk emitting the big blue bump (BBB) can be written as 



(Collin 2001), where T§ is the temperature of the disk in units of 10 5 K, M§ is the mass 
of the black-hole in units of 10 8 M©, R is the radial distance from the SMBH, Rg is the 




(2) 
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Fig. 16. The observed photon index in the 2 — 10 keV band, derived from the time 
resolved spectroscopy of 2000 (filled circles) and 1994 (open circles) observations, plotted 
against the observed count rate in the same energy band. The count rate is the average 
of the GIS2 and GIS3 count rates. Also shown is the best-fit function of the form Tx = 
T - KF- S . The best-fit parameters are T = 3.24 ± 0.97, S = 0.39 ± 0.32. 



Schwarzschild radius, Lx and LEdd are the X-ray and Eddington luminosity, respectively. 
For NLS1 galaxies, T ~ 10 6 K, M ~ 10 6 M©, thus n e ~3x 10 12 cm" 3 for R = WR S . 

If the disk thickness H ^ -± 10 12 cm ~ 3i?s, then the disk is optically thick for 

Thomson scattering. Thus for super Eddington rates, the disk can be optically thick 
and the high energy photons from the corona incident onto the disk will lose energy by 
direct Compton scattering and the disk electrons will gain energy. Due to the Coulomb 
interaction, electrons will quickly thcrmalizc thus increasing the disk temperature and 
hence increasing the soft hump emission. The soft hump emission of NLS1 galaxies, 
therefore, may consist of intrinsic disk emission, reprocessing of hard X-ray emission by 
photoelectric absorption and by direct Compton scattering. 
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8.3. Correlation between the power-law flux and the soft hump emission 

The soft-band flux and the power-law flux change either simultaneously or with short 
time delay (see Fig. ||). Our time resolved spectroscopy has revealed that the soft hump 
flux in the 0.7 — 2.0 keV band changes by larger factor (~ 4.7) than the change (by a 
factor of ~ 3.4) in the power-law flux in the 2 — 10 keV band on a timescale of ~ a week. 
However, on a timescale of 20000 s the power-law flux changes by a larger factor (~ 6) 
than the change in the soft hump flux (by a factor of ~ 2.6) (See Fig. |^ and Sect. ||). 
The above trend suggests that it is the power-law component that is responsible for the 
most rapid (^ 1000 s ) variability while the soft hump dominates the longer timescale 
(<; a week) variability. Thus the 500 s variability reported here and 800 s variability of 
IRAS 13224-3809 reported in Boiler et al. (1993) in the ROSAT band could be entirely 
due to changes in the power-law component. However, better signal-to- noise data over 
a broad energy band, for example with XMM-Newton is required to varify the above 
idea. Boiler et al. (1993) rejected the standard thin accretion disk model, in spite of 
the good-fit to the ROSAT PSPC data, on the ground that the standard disk emission 
cannot produce the observed soft X-ray variability as the shortest timescales (e.g. thermal 
timescale, sonic timescale) possible for standard thin disks are longer than the observed 
variability timescale by a factor of 2 — 3. This apparent problem can be resolved if the 
observed rapid variability in the ROSAT band is due to the changes in the power-law 
component alone. If the heating of the corona is by magnetic reconnection (Merloni 
& Fabian 2001), the variability timescale could be as short as the coronal dissipation 
timescale which is given by t(R) ~ i <smq where R is the size of an active region 
(see Merloni & Fabian 2001) and u<n a is the dissipation velocity, a few. Thus 

t(R) il a few x 100 s and rapid variability of the power-law flux can be produced via 
inverse Compton scattering of disk photons in the active regions. An increased flux of the 
power-law component would further heat the disk due to increased irradiation. However, 
detailed time dependent accretion disk-corona models are required in order to understand 
the correlation between the soft hump emission and power-law flux and the variability 
amplitudes at different timescales. 



9. Summary 

1. On a 10-day baseline, the 0.7 - 10 keV band flux of IRAS 13224-3809 shows trough- 
to-peak variation by a factor > 37 when sampled using 500 s bins. The hard-band 
(1.3 — 10 keV) and soft-band (0.7 — 1.3 keV) fluxes present trough-to-peak variations 
by factors > 25 and ~ 20, respectively when sampled using 5000 s bins. 

2. The intensities in the soft and hard bands are strongly correlated and are due to 
changes in both the power-law flux and soft hump emission simultaneously or perhaps 
with short time (~ a few x 100 s ) delay. 
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3. The fractional variability amplitude is variable in the soft and total bands but not in 
the hard band. None of the variability amplitude is correlated with the X-ray flux. 

4. We confirm the presence of a "soft hump" above the power law at energies below 
< 2 keV. The power-law component contributes only ~ 32% of the flux to the soft 
hump in the 0.7 — 1.3 keV band. The soft hump component shows flux variations 
down to timescalcs of 1 day, ranging by a factor of 4.7 in the 0.7 — 2 keV band. The 
soft hump flux appears to be correlated with the power-law flux. 

5. The mean photon index is 2.11^Q gg. Time resolved spectroscopy reveals significant 
changes in Tx, AI"x = 0.57 ± 0.21 on a timescale of 5 days. However, day-to-day 
variations are not significant. Variations in the power-law flux on times of ~1 day are 
not due to changes in the power-law slope. 

6. In our daily sampling, the photon index seems to be correlated with the flux of the 
soft hump, but due to poor signal-to-noise of the data, a firm conclusion cannot be 
made. 

7. Flux-selected spectral fits reveal that at higher flux levels the power law becomes 
steeper and the soft hump and power-law flux higher. The photon index changes by 
AFx ~ 0.4 while the soft hump flux changes by a factor of 5.6 and the power-law flux 
by a factor of 3.2 between the low and high flux states. 

8. The 1994 ASCA observation also reveals increase in the Tx by > 0.2, the soft hump 
flux by a factor of ~ 6 and the power-law flux by a factor of <~ 2.8 from the low to 
high flux states. 

9. An Fe Ka line is not detected from IRAS 13224-3809. The 90% upper limit to the 
equivalent width of Fe Ka is 400 eV. 
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